Semiclassical limits to the linewidth of an atom laser 
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We investigate the linewidth of a quasi-continuous atom laser within a semiclassical framework. 
In the high flux regime, the lasing mode can exhibit a immber of undesirable features such as density 
fluctuations. We show that the output therefore has a complicated structure that can be somewhat 
simplified using Raman outcoupling methods and energy-momentum selection rules. In the weak 
outcoupling limit, we find that the linewidth of an atom laser is instantaneously Fourier limited, 
but, due to the energy 'chirp' associated with the draining of a condensate, the long-term linewidth 
of an atom laser is equivalent to the chemical potential of the condensate source. We show that 
correctly sweeping the outcoupling frequency can recover the Fourier-limited linewidth. 
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I. INTRODUCTION 

Optical lasers have found broad application in preci- 
sion measurements that address questions both funda- 
mental and applied in nature. In many cases, we expect 
to be able to perform such experiments more effectively 
and to higher precision with atom interferometry [l], 
Ultra-cold ensembles of thermal atoms are already uti- 
lized in interferometric systems where they have been 
demonstrated to compete with the best inertial and grav- 
itational measurement apparatus available P, H, 0) B ■ 
For this reason, there is significant interest in the produc- 
tion of a coherent source of ultra-cold atoms from a Bose- 
Einstein condensate (BEC) for applications in precision 
measurement and metrology 0, Q ■ A free-space BEC 
atom interferometer in the Mach-Zchnder configuration 
has been demonstrated to produce 100% contrast at the 
output port 0, and there are proposals to put atom in- 
terferometers into space [l^]. Due to trap stability and 
mean-field effects, precision experiments will most likely 
need to be performed with the low density, untrapped 
atomic beam rather than in more dense atomic sources 
such as full Bose Einstein condensates • There is also 
interest in producing and measuring non-classical quan- 
tum states of atomic beams [ll|, [13, [l^l- All of these 
proposed applications will require a spatially stable atom 
laser beam with good first order coherence. 

In this paper we investigate the key spatial properties 
of a quasi-continuous atom laser. In the cases where the 
outcoupled atom laser is stable, we focus on the linewidth 
of the output spectrum as the key measure of the first- 
order coherence of the beam. We begin by highlight- 
ing the particular importance of linewidth for dispersive 
fields such as atoms. In Sec. IIIII we introduce our model, 
and we calculate the properties of the output in various 
limits in Sec. IIVI In Sec. W\ we describe methods for 
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reducing the atom laser linewidth. 



II. THE IMPORTANCE OF LINEWIDTH 

In a precision interferometric measurement made at 
the shot noise limit, all that is theoretically required of 
the wave source, whether it be a source of matter waves 
or light, is that it have high flux. In principle, classi- 
cal source fluctuations in frequency and phase can be 
removed through good interferometer design. A long co- 
herence length, equivalent to a spectrally narrow source, 
is not required if the path length difference in the inter- 
ferometer is less than the coherence length. In principle, 
mode matching on the output beam splitter of an inter- 
ferometer can be performed as well on a complicated spa- 
tial mode as a simple one and a highly divergent beam 
can be collimated with lenses. In practice, however, if 
an interferometer is to operate at the shot noise limit, 
none of this is true. The shot noise limit for a high flux 
source is difficult to achieve, and it is essential to have a 
spectrally narrow, classically quiet, low divergence beam 
with the minimum transverse structure in both phase 
and amplitude. For these reasons, typical precision op- 
tical measurements use classically quiet lasers operating 
on the TEMoo mode. 

These points apply equally well to optical and atom 
lasers, suggesting that in the absence of a perfect experi- 
ment, a spectrally pure output beam is highly desirable. 
The difference is that unlike in an optical interferometer, 
the linewidth of a matter wave interferometer can still 
be critical even for a perfect experiment. To see this we 
consider the case of an equal path length Mach-Zehnder 
interferometer and an atomic beam that has an average 
momentum hk and a momentum spread of h 5k. We now 
introduce some disturbance in one of the arms, corre- 
sponding to the effect we wish to measure, and model it 
as a step function potential of width L and height Vq. 
This is shown schematically in Figure [1] 

Assuming that 5k <^ k, it is easily shown that due to 
the dispersive nature of the atomic beam a phase differ- 
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FIG. 1: Two matter waves with wave vectors k and k + 5k 
are incident on a step potential of height Vo- 
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builds up between the k and k + Sk components of the 
incident beam. The phase difference Acj) represents an 
inherent uncertainty in the phase resolution of the inter- 
ferometer, regardless of what measurement technique is 
used at the output ports. This uncertainty arises purely 
from the finite linewidth of the atomic beam, and will 
exist no matter how accurately the path lengths of the 
interferometer arms are matched. 

Now consider the case for an optical interferometer, 
again using spectrally broad beams, and again with path 
lengths perfectly matched. If a potential is introduced 
into one of the arms, it corresponds to a change in the 
refractive index over that region. Provided this change 
is not dispersive, all the components in the beam will see 
the potential simply as an increase in path length — cru- 
cially, the same increase. Consequently the interferome- 
ter can be nulled and brought back to correct operation 
by physically readjusting the distance in one of the arms. 
As indicated in Eq. ([T]), this readjustment is not possible 
in an atom interferometer. 

In general, generating an atomic source with low 
linewidth also requires the beam to have well-controlled 
spatial properties in other respects, which makes it a 
good choice of metric for the first order coherence of 
atomic sources. Applications that do not specifically 
benefit from the narrow linewidth of an atom laser will 
still tend to benefit from the associated controlled spatial 
mode. 



atoms into the untrapped state where they no longer ex- 
perience the confining potential and are free to leave the 
trap. 

There are two common outcoupling methods, both of 
which require that the atoms of the BEC are in a spe- 
cific angular momentum substate, for example mp — 1, 
meaning they can be confined by a magnetic trap. The 
two methods are shown schematically in Figure [21 

In the first outcoupling scheme, an external rf field is 
applied that flips the atoms into an mp =0 state that 
does not see the trapping magnetic field. Consequently 
the atoms fall from the trap under the influence of grav- 
ity, creating a semi-directed beam of coherent atoms: 
the atom laser. In the second scheme a Raman outcou- 
pling method is used, where two optical fields transfer the 
state of the atoms from trapped {mp = 1) to untrapped 
{mp — 0) via a third intermediate level. The use of two 
optical fields means that a significant momentum kick 
can be imparted to the untrapped atoms as they leave 
the condensate, leading to an atom laser with superior 
properties such as a higher flux and higher brightness 11411 
and the ability to give the beam directionality [13, Il5j |. 
In addition, the Raman scheme allows for the possibility 
of creating non-classical states of the beam [ij, [l^ . 
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FIG. 2: An atom laser based on a) Raman outcoupling and 
b) rf outcoupling. In both cases trapped atoms in state |1) 
are transferred to an untrapped state |2) via electromagnetic 
fields. The Qij represent the Rabi frequency of the applied 
fields and the Aij represent detuning from resonance. 



III. MODEL 

In the most general terms, an atom laser requires cou- 
pling atoms out of a BEC into a coherent beam. The 
most common way to accomplish this is to use a state- 
changing outcoupling method, where the atomic species 
making up the BEC has at least two separate internal 
states — a trapped state and an untrapped state. The 
BEC consists of atoms in the trapped state, which feel 
some confining potential that keeps them localized and 
Bose condensed. Some external perturbation is then ap- 
plied to the BEC which flips a portion of the trapped 



Regardless of whether an rf or a Raman outcoupling 
scheme is used, the second-quantized Hamiltonian de- 
scribing the system can be written most generally as 

-ff = y (^Htiup + Hhcam + Hint^ (fir (2) 

where -&trap describes the atoms in the trap, ffboam de- 
scribes the atoms in the atom laser beam, and iJjnt de- 
scribes the outcoupling process. For the purposes of this 
paper we assume an isotropic harmonic trapping poten- 
tial for the trapped atoms, which enables us to write the 
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terms of Eq. ^ with greater specificity as 



IV. LINEWIDTH CALCULATIONS 
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i^beam = *I (^"^V^ - M + j ^„ (4) 

Hint = -?i(ri(r)*,t^'t + r!*(r)*|4'„ 

+[/t„^I'l^t,I,„,I,,) (5) 

where ^'t(r) and *Pu(r) describe the trapped and un- 
trapped matter fields respectively, ujt is the harmonic 
trapping frequency, and the Uij are nonlinear potentials 
arising from atom-atom collisions. For rf outcoupling the 
detuning is given by 6 = Art and for Raman outcoupling 
it is given by 
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The specifics of the outcouphng strength 0(r) depend 
on the outcoupling method. In the rf outcoupling case 
n = riff and has no position dependence. With Ra- 
man outcoupling, however, is position dependent and 
is given by 



Ai3 



(7) 



where kp = k2 — ki is the momentum kick imparted to 
the outcoupled atoms from the two optical beams. 

We assume that the atomic gas is sufficiently cold and 
dilute, so only binary collisions are relevant and the non- 
linear potentials are defined by 



Uij — 47r ?i^ay/m. 



(8) 



where aij is the s-wave scattering length between atoms 
in state |i) and state 

When the quantum statistics have no effect on the dy- 
namics of the mean field, we can use the Gross-Pitaevskii 
(GP) equation to describe the atom laser. The coupled 
CP equations arising from Eqs. ([2]) - ([5]) are given by 
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dt 
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+Utu\^tn^Ju-nn*{r)^Pt. (10) 



These equations can be solved numerically in one, two or 
even three dimensions, depending on the spatial resolu- 
tion required and computational resources available. 

As we have included position dependence in the matter 
fields, the effective Hamiltonian describes the full mul- 
timode nature of the problem, and also includes non- 
Markovian effects. 



Unlike the optical laser, atomic fields do not have a 
simple proportional relationship between the energy and 
momentum spectra, and hence the linewidth in one does 
not translate trivially to a linewidth in the other. In 
free space, the distinction is largely irrelevant as both 
spectra are static, but this is not always true. In most 
current experimental atom lasers, atoms are outcoupled 
from a trap and allowed to fall under gravity. As they 
fall, they gain kinetic energy at the expense of potential 
energy. This results in a kinematic compression effect, 
leading to a narrower spread in momentum the further 
the atoms fall, while the energy spread remains constant. 
The energy spread of the beam is thus a more stable mea- 
surement of the beam's linewidth than the momentum 
spread. 

Before numerically solving the full equations ^ and 
(llOp , it is worthwhile studying a number of simplified ver- 
sions of the problem in order to extract as much analytic 
insight as possible from the problem. To this end, we will 
first examine the problem in the case where the atomic 
nonlinearities in the condensate are negligible, and as- 
sume the condensate remains single mode. We then ex- 
amine the case where the condensate is allowed to be 
multimode, but still linear. Finally we consider the fully 
general case numerically. 



A. Single-mode condensate with no nonlinearity 

To begin, we assume the condensate begins in a single 
mode and remains single mode due to the outcoupling be- 
ing weak enough such that there is negligible back action 
on the condensate. This is the simplest possible case, and 
will result in the minimum possible linewidth achievable 
in a non-pumped atom laser. 

We will also assume that the condensate has no nonlin- 
ear interactions, which means we set Uu = Uuu = Utu = 
in Eqs. ([9]) and pO)) . This represents a regime where 
the condensate contains few atoms, the trapping poten- 
tial is weak, or the s-wave scattering length of the atomic 
species is small. 

In this situation one might naively assume that the 
linewidth of the atom laser is the power-broadened 
linewidth of the atomic transition. To see this, one can 
consider the outcoupling process as a classical picture of 
"atoms rolling down a hill" , where the untrapped atoms 
are created on a potential hill arising from the fact that 
the potentials the trapped and untrapped atoms see are 
different. In this naive picture the atoms appear on the 
potential hill with a spread in position given by set- 
ting the energy spread equal to the power broadened 
linewidth, as depicted in Figure O This picture is false, 
however, as the wave-like nature of the atoms can cause 
destructive interference of some energies and constructive 
interference of other energies. It is necessary to take the 
wave-like nature of the atoms into account to accurately 
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describe the energy spread. 



^(co + 6) - 
h (0 



\ 

J 


/ 

BEC / 










— ^^^^ 



Z (arbitrary units) 

FIG. 3: Naive description of tlie energy spread of outcoupled 
atoms due to a "balls on a hill" model. The atomic transition 
has a resonant width 5 such that the coupling happens over a 
region Az, corresponding to a change in potential of h,5. This 
leads to a spread in energies of Ti5 for the outcoupled atoms. 
This model neglects the wave-like nature of the atoms, and 
that there is interference between the different energies. 



The Hamiltonian given by Eq. ([2]) consists of terms de- 
scribing the trapped matter field, the untrapped matter 
field, and the coupling between the two. We denote the 
ground state energy eigenfunction of the trapping Hamil- 
tonian i?trap as (j^tix) and the energy eigenfunctions of 
the beam Hamiltonian i/bcam as where q is any 

convenient continuous parameter that can label the en- 
ergies of the atom laser beam. The eigenvalues of (/>t(a;) 
and 4>u{q,x) are Uloq = hLOt/2 and fiLj{q) respectively. 

We can now expand the wave functions of the trapped 
and untrapped matter fields as 



i^uix^t) = / P{q,t)4>uiq,x)dq. 



(11) 
(12) 

I and 



icxQ = UJoaQ-ft / A{q)P{q,t)dq (13) 

J — OO 

t${q,t) - iu;iq)-5)P{q,t)~n*A*{q)ao{t), (14) 



Under the approximations described above, Eqs. 
(ITOl) become 



where A{q) — J_^(f>f{x)<j)u{q,x)A{x) dx, and A{x) — 
il{x)/n represents the spatially-dependent part of the 
electromagnetic field coupling the atoms out of the trap. 



and in 



In the case of Raman outcoupling A(a;) 
the case of rf outcoupling A(a;) = 1. 

Eqs. lfT3|) and can easily be solved numerically, but 
we first derive an approximate analytic solution to gain 
insight into how the lincwidth scales with various param- 
eters. 



By making the transformation 



Piq,t) = /3(g,i)e^(-(')-^'*, 



(15) 
(16) 



we obtain 
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A{q)p{q,t)e"^'^^l'>'dq (17) 



iP{q,t) - -l]*A*(g)5o(i)e"'^"^'^*, 



(18) 



where Aa;(q) — loq — {uj{q) — 5). Formally integrating Eq. 
(ITSl) and assuming the initial state of the output field is 
vacuum we obtain 

P{q,t)^in*A*{q) f ao{t')e"^'^'^'^^'' dt' . (19) 

Substituting this result into Eq. ^T7\ gives 

ao = [ f \A{q)\^a,{f)e^^-'-'^-^'-'-''^^du:dt'. 

(20) 

To proceed we make use of the fact that in the weak 
outcoupling regime the momentum spread of the output 
is much narrower than the momentum spread of the con- 
densate. This means that due to momentum conservation 
we are only selecting atoms with a narrow range of mo- 
menta from the condensate, and thus can assume that 
the form of A{q) is flat over this range, allowing us to 
replace A{q) with A{qQ) and take it outside the integral. 

Similarly, due to energy conservation, the energy of 
the outcoupled atoms will be centered around Ti{uj{) -I- 5), 
meaning most of the dynamics will occur at frequencies 
close to + Provided dq/duj(q) is slowly varying close 
to ujQ + 6, it is a valid approximation to replace dq/duj 
with dq/duj\uja+s- For free space dq/doj oc which 
means that for a Raman transition with a large momen- 
tum kick it will be approximately constant. In the case 
of a gravitational potential, dq/du is constant, so the 
approximation is exact. Using these approximations we 
obtain 
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m'\Mqo)\' ^ 
duj 
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aoit')e''^^^''^^'~''^dwdt' 



Juj— — oo 



where 



j = 2nmA{qo)\'^ 
duj 



(21) 



(22) 



As ao{t = 0) = VNq, where A^o is the number of atoms 
in the condensate at i = 0, the solution to Eq. (pij) is 



<5o(0 = vAoe" 



(23) 
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Consequently, the condensate number N{t) = \ao{t)\'^ 
will decay exponentially according to N{t) = Nqb'^*. 

The spectrum of the output (i.e. number of atoms per 
mode q) can now be obtained from Eq. (fTO|) . We find 



\f3{q,t)\^ = mA{q)\^N,F{Aio,t) 



where 

F{ALj,t) 



1 - 2cos(Aa;i)e-i* +e-'^*' 



(24) 



(25) 



This spectrum is plotted in Fig. 3) 
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FIG. 4: F{Auj,t) for 7 = 1 Hz. As i ^ oo, F{Auj) asymptot- 
ically approaches a Lorentzian of FWHM 7. 



In the weak coupling limit, when A{q) varies with uj 
much more slowly than F{Au!,t), as t 00 the spec- 
tral power density of the output beam |/3((j')P becomes a 
Lorentzian with a full- width- at-half-maximum given by 
7. Thus, in the long time limit, the fundamental limit 
to the linewidth of an atom laser is related to the time 
it takes to drain the condensate: 7 = tJ^^^„, where 7 is 
the spectral linewidth (measured in rad s~^), and Tdrain 
is the 1/e drain time of the condensate, measured in sec- 
onds. This result says that the linewidth of a pulse of 
atoms coupled out of the condensate sufficiently weakly is 
only limited by Fourier arguments, although it should be 
noted that most atom laser experiments are in a stronger 
coupling regime. The weak coupling limit gives us the 
same separation of timescales found in optical cavities, 
and thus the spectrum is identical to the spectrum of 
photons draining out of an optical cavity. As the average 
flux from our atom laser is Tav = Nq/ At, (where At is 
either the drain time of the condensate, or an artificially 
imposed cut off time of our atom laser pulse), and our 
spectral linewidth is always limited to Soj > 1/At, we 
obtain the inequality 
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TABLE I: Comparison of condensate drain time 1/r and the 
long-time linewidth of an atom laser for different values of 
fl, and ko- 7n and 7a represent the linewidths for the ana- 
lytic theory and numerical simulation respectively. The table 
shows close agreement between our approximate analytic re- 
sult and our numeric calculation. The entries marked with 
an asterisk display poor agreement between the analytic and 
numeric results. This is because the coupling is sufficiently 
large such that the approximation made in Eq. (|2H) is invalid. 



relating the average flux and spectral linewidth of an un- 
pumped atom laser. 

As a number of approximations were used to obtain 
Eqs. dMl) and we also solved Eqs. ^ and dH]) nu- 
merically and compared the results. For our numerical 
model we took parameters typical to a ^^Rb atom laser 
such as one described in [1J|, and chose ko = 10'' m~^, 
LUt = 50rads~-^, and O = 50rads~^. We assumed out- 
coupling into free space, and with the free space disper- 
sion relation Eq. (22) becomes 



7 



m 1 



(27) 



< No 



(26) 



TableUshows the comparison between the analytic the- 
ory and the numerical simulations. Overall, there is good 
agreement, allowing us to use Eqs. ([M)) and ([77|) with 
some confidence. 

The key point of this semiclassical analysis is that 
for an ideal, single-mode, non-pumped atom laser, the 
linewidth is given by the inverse of the drain time of the 
condensate. Consequently it can be made as narrow as 
desired by reducing the outcoupling strength arbitrarily. 
The trade off is that this arbitrarily narrow linewidth 
comes at the expense of reduced flux. 



B. Multimode condensate with no nonlinearity 

We now consider a more realistic model, where the 
condensate is not constrained to remain in a single mode, 
although we still assume the nonlinearities are negligible. 

We proceed as in the previous section, except now we 
allow the condensate to be multimode, and for simplicity 
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treat both the condensate modes and atom laser modes 
as discrete. The wave functions of the trapped and un- 
trapped atoms are now expanded as 



ipt{x,t) = y^Q„(^)(/)f„(x) 

n 

1puix,t) = ^ (3n{t)(f>un{x) 



(28) 
(29) 



where (j)tn{x) and (j^unix) arc the nth eigenstates of the 
Hamiltonians iJtrap and iJbcam respectively. We denote 
their eigenvalues by hujtn and fiUun- Proceeding as be- 
fore, we find that the equations of motion in an appro- 
priate rotating frame are given by 



-i{uJtn — {^^um—S))t 



(30) 



(31) 



where A„™ = (/ij*„(x)0„„(x)A(a;) da;. 

When the coupling is weak, the phase rotation of 
am{t) and Prn{t) is approximately at zero frequency. This 
means that over long times the only significant contribu- 
tion to the growth of /3m (i) is from the trapped mode a„ 
with frequency 



(32) 



as all other modes will on average cause no net growth 
on time-scales much larger than r = (tutn — i^mn ~ S))^^- 

In the case of rf outcoupling, the momentum kick to the 
outcoupled atoms is negligible, so the atoms retain the 
momentum they had when they were in the trap. This 
means the energy conservation relation ((32|) is the only 
condition that must be satisfied when considering the 
output spectrum, and consequently the energy spread of 
the output is now related to the energy spread of the 
condensate via the magnitude of the matrix elements 
Anm- Thus the spectrum of the atom laser will essentially 
mimic the spectrum of the condensate, with the energy 
peaks in the output beam at frequencies Wt„ -I- S corre- 
sponding to different energies in the condensate, moder- 
ated by the magnitude of A„m, with each peak broad- 
ened such that it is the Fourier limit of the outcoupling 
time, as discussed in the previous section. It is therefore 
clear that in the rf case any dynamic fluctuations in the 
BEC will result in an atom laser with a broader spectral 
linewidth than if the BEC were single mode. 

The situation is more complicated in the case of a Ra- 
man outcoupling scheme with a large momentum kick. If 
we assume that our condensate is initially in the super- 
position 



'ipt{x) = ^a„0t„(a;). 



(33) 



then, for long times, only the states that satisfy the en- 
ergy resonance totn = ^uj — ^ will be present in the out- 
put, just as in the rf case. However, now there is an 



additional constraint arising from the dispersive nature 
of the atoms. Taking free space as an example, we have 



2m' 



(34) 



huj (x) . Thus the output 
will only contain momentum states 



where hkj is the momentum of < 



kj — 



2m{5 + LOn) 



(35) 



However, conservation of momentum demands kj = k + 
ko, where k is the initial momentum of the atoms. The 
source of outcoupled atoms will then be the component 
of the condensate with momentum 



kn 



2m{5 



kn. 



(36) 



The relative amplitude of the component of the output 
with momentum kj is then proportional to the value of 
An{k), the fc-space representation of (j)t^n{x)^ evaluated 
at kn- This can lead to interesting effects. As an exam- 
ple, we choose our trapped system as a harmonic oscil- 
lator. When fco is small, the spacing between the output 
k modes is of order Afcj ~ ^^.mujtl'h^ which is about 
the spacing of the 'lobes' in the k space representation of 
^tn{x)- When /eq is large, Afc^ ~ 0, so the outcouphng al- 
ways happens close to the center of the momentum space 
wave function. As An(k) are Hermite Gaussians in the 
case of the harmonic oscillator, An(k 0) = for n 
odd, so there is no outcoupling from odd modes. This is 
represented in Fig.O 

The relative intensity of each k component in the 
output is proportional to |Q;„p|yl„(A:„)p. Fig. [6] shows 
|^n(^n)P for different values of k^. Fig. [5] was checked 
against a multimode Gross-Pitaevskii simulation for the 
first five modes. Close agreement was found for the rela- 
tive amplitudes, although the momentum resolution was 
insufficient to accurately resolve the difference in mo- 
mentum for each peak. For small fcg, all modes are 
present in the output, with the relative intensity of each 
mode decreasing due to the spreading out of Anik). For 
large fco, only even modes are present in the output. 
This is due to fc„ falling at nodes of An(k) for the odd 
modes. In the cross-over regime (when fco is of order 
^JrrujJtjK) , there is complicated structure in the output 
with some even and some odd modes severely attenu- 
ated in the output. For *^Rb, at a typical trapping fre- 
quency (wt = 50rads^^), this cross over occurs at around 
fco ~ 4 X 10^ m~^, which is much less than the maximum 
recoil of fc « 1.6 x 10^ m~^ achievable with a two-photon 
transition, using light of wavelength A = 780 nm. This 
suggests that this effect should be observable in exper- 
iments, although our theory has neglected the atomic 
interactions, which will complicate the effect. 



7 



L=10 m 



n=0 



«=3 



/t„=10 m 



(1=0 




«=3 



11=8 



/ \ 



L=10 m 



n=0 



(1=3 



y I 



(1=8 



-2-1 ' 2-2-1 1 2-2-1 1 2 



k-k^ (I'm-'] 



1 r 

0.5 



(a) 



° o o 



o o oooooooooo 



^ 2 4 6 

I 1, 

0.5 

CD ° ° o ° o 

> qI ^ ° i ° i ~ ° 1 o ' 

01234567 
P 1 



10 12 14 16 18 20 



X 10 



(b) 



9 10 5 
xlO 



(c) 



0.5 







10 



15 



k-k^ (m') 



xlO 



FIG. 5: |A„(fc)|^ for n = (top), n = 3 (middle) and n = 8 
(bottom), kn — 1/(5 + ujri)2m/'h — fco, the place in the k- 
space wave function from where resonant outcoupling occurs, 
is indicated by a vertical black bar in each case. For fco = 
10^ m~^, kn follows the largest 'lobe' of An{k) for increasingly 
excited states. For fco = 10* m~^, fc„ remains approximately 
in the center of An ik) . 5 = Tik^ /2m — ujto was chosen such 
that the outcoupling was perfectly on resonance for the zero 
momentum component of the ground state. 



FIG. 6: The relative intensities of the momentum components 
of the atom laser beam corresponding to the first 20 conden- 



sate eigenmodes for (a) fco 
-1 



lO"" m" 



(b) fco = 10" m"^ and 
(c) fco = lO'^m'V For small fco (fco = lO'^m"^), the relative 
intensity of each mode decreases due to the spreading out of 
An{k). For large fco (fco — 10* m~^), only even modes are 
present in the output. This is due to kn falling at nodes of 
An{k) for the odd modes. In the intermediate case, there is 
complicated structure, with some odd modes and some even 
modes being attenuated. 



C. Multimode condensate with nonlinear 
interactions 

We now turn to an analysis of the complete problem 
and allow the nonlinear interactions in the condensate to 
be significant, which is the case in many experimentally 
realizable atom lasers. 

The existence of nonlincarities makes a difference to 
a number of properties of the condensate and the atom 
laser, affecting things such as mode shapes, memory func- 
tions and classical density fluctuations. The change that 
is most relevant to the linewidth, however, is the fact 
that the energy of the condensate is now dependent on 
the number of atoms in the condensate. 

To understand this, we note that when an atom in the 
condensate is flipped from a trapped to an untrapped 
state, it experiences a mean field potential that depends 
on the density distribution of the condensate. As the 
density of the untrapped field is much less than 

that of the trapped field inside the condensate, it is clear 
from Eq. (jlOp that this mean field potential is given by 

V^i^Uuu\M'r)\\ (37) 

where Uuu is defined by Eq. ([8]) . As density is always pos- 
itive and increases towards the center of the condensate, 
this results in a repulsive force on the untrapped atoms 
that accelerates them out of the BEC. 

As an example, we consider the Thomas-Fermi limit. 



where the nonlinear energy of the condensate is consider- 
ably larger than the kinetic energy and the wave function 
for the condensate can be found analytically. Assum- 
ing a harmonic trapping potential with frequency cJt, the 
atomic density in the condensate is given by 

= ^ [^^{N) - \muy - mgz^ (38) 

where 

Eqs. ([57)1 and (|38p show that after the atoms have been 
flipped into an untrapped state, they slide down a 
quadratic potential hill of height fi, giving them a kinetic 
energy as they leave the condensate. 

Assuming the condensate is not pumped, its atom 
number will inevitably reduce during the outcoupling 
process. As the number of atoms in the condensate 
falls, the the chemical potential also falls, resulting in 
a situation where atoms outcoupled later in time will 
have a smaller kinetic energy as they leave the conden- 
sate compared to atoms outcoupled at an earlier time. 
Consequently, the atom laser beam will consist of atoms 
with a wide spread of energies; a spread that can be as 
large as n{No) if all the atoms are outcoupled. 
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To demonstrate this effect we numerically solve the full 
Eqs. © and PH)) for a situation with experimentally real- 
istic parameters and a large nonlinearity. The equations 
were solved in one dimension only, using a dimensional 
reduction procedure where the nonlinear potentials Uij 
were scaled by a transverse area corresponding to the 
cross sectional area of the beam [l^ • This not only makes 
the computation far more tractable, but also removes ad- 
ditional complications that obscure the linewidth such as 
part of the mean-field kick being transferred into trans- 
verse modes of the laser. The simulation couples into free 
space rather than a gravitational potential, meaning the 
momentum space wave function of the beam can be used 
to give the linewidth. 

Figure [7] shows the result of the simulation, displaying 
snapshots of the beam's momentum space wave function 
at various points in time. Initially the momentum is cen- 
tered at fccont = y^/jQ + 2m/i(A*'o)/?i^, and the linewidth 

(i.e. momentum spread of the beam) begins to narrow 
in accordance with the Fourier argument laid out in Sec- 
tion HVAl After enough atoms have been outcoupled to 
significantly change the chemical potential, however, new 
atoms appearing in the beam have lower and lower en- 
ergies, resulting in the line center being "chirped" . Con- 
sequently the effective linewidth becomes ever broader, 
ultimately spanning all momentum states between fco and 
^cent when the entire BEC has been drained. 

The complicated structure seen in Figure [7] arises from 
the fact that an atom laser beam is a complex field, and 
as the frequency of condensate phase evolution is chang- 
ing, there can be destructive interference between atoms 
emitted with a particular energy and atoms emitted at a 
later time with the same energy but out of phase. 



V. METHODS TO REDUCE LINEWIDTH 

In this section we review possible methods for reducing 
the linewidth from an unpumped atom laser in various 
parameter regimes. 



A. Weak outcoupling 

The most obvious way to minimize the linewidth of 
an atom laser is to outcouple extremely weakly, as this 
increases the drain time, effectively without limit. Thus, 
provided the condensate nonlinearities are negligible, the 
Fourier arguments in Section IIV Al demonstrate that the 
linewidth can be made arbitrarily narrow. 

In the case where condensate nonlinearities are not 
negligible, weak outcoupling still succeeds in reducing the 
linewidth as the condensate will undergo almost no deple- 
tion, meaning the chemical potential is static and there 
is no chirp of the line center. However, unlike the case 
where nonlinearities can be ignored, we cannot weakly 
outcouple all the atoms — we must ensure that over the 
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FIG. 7: Momentum space density of an atom laser beam after 
a) 20 ms, b) 100 ms, c) 400 ms, and d) 1000 ms of outcoupling. 
After 1000ms 65% of the atoms in the condensate have been 



outcoupled. 

m, 



Parameters: A^o 
fco = 3.2 X 10*^ m" 



10" 



150 rads 



entire duration of the experiment the change in chemi- 
cal potential is less than the minimum linewidth we are 
willing to accept. In the case of a strongly nonlinear 
condensate, if we require the temporal linewidth to be 
less than Siv, then by Eq. ([55)1 we must ensure that the 
number of atoms outcoupled from the condensate is less 
than 



AN < 



mujfh 



2/5 



(40) 



where ujt is the harmonic trapping frequency and A^o is 
the number of atoms initially in the condensate. 

This approach can be arbitrarily effective if high flux 
is not important. We solved Eqs. (O and (fTU)) for a 
highly nonlinear system, and examined the linewidth 
of the beam over time. Over the simulation approxi- 
mately 10 atoms were removed from the condensate, cor- 
responding to extremely weak outcoupling. The results 
are shown in Figure[Hl In the long time limit the curve 
shown in Figure [8] is linear with a slope of -1, indicating 
that linewidth is inversely proportional to the outcou- 
pling time. This agrees with the Fourier arguments our 
approximate single-mode linear theory predicts. 

The difficulty with weak outcoupling is that high flux is 
one of the more desirable qualities in a laser. What have 
demonstrated here is that in the ultra-low flux limit the 
atom laser can be regarded in some sense as having a 
very narrow linewidth with a slowly moving line center. 
In practice, this flux limit will make any experiments im- 
practical, so we now consider a method to achieve narrow 
linewidth without sacrificing flux. 
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FIG. 8: Linewidth narrowing as a function of outcoupling 
time. Parameters: A'^ — 10^, ujt = 250rads~^, a = 3xl0~^m, 
fco = lO'^m^V 



B. Chirp compensation 

The source of the drift of the hne center is the mean 
field potential that untrapped atoms experience as they 
leave the condensate. If we ignore gravity, this mean field 
potential can be found from Eqs. ((37)) and (|38p and is 
given by 

Vmf = Kt) - Imcjy, (41) 

where we have taken Uu — Uuu and allowed the chemical 
potential time dependence to take into account conden- 
sate depletion. If we choose the outcoupling point to 
be the center of the condensate, atoms acquire the full 
chemical potential ^{t) worth of energy on their way out 
of the condensate. However, as the condensate depletes, 
pL{t) decreases, meaning atoms outcoupled later have less 
energy, broadening the linewidth. 

The solution is to begin outcoupling from a point away 
from the center of the condensate, so that initially atoms 
do not acquire the full worth of energy as they slide 
down the potential hill. If we then move the outcoupling 
point back in towards the center of the condensate as 
it depletes, it is possible to ensure that atoms outcou- 
pled later acquire the same amount of energy as those 
outcoupled earlier, thus removing the chirp effect. 

This shifting of the outcoupling point can be accom- 
plished by making the two-photon detuning 5 time de- 
pendent. If we wish to initially begin outcoupling atoms 
from a distance tq from the minimum of the magnetic 
trap (which will coincide with the center of the conden- 
sate if gravity is ignored), and sweep this point towards 
the center of the trap in such a way that atoms always 
leave the condensate with the same energy, then, since 
the energy of atoms as they leave the condensate is 

£^out = S{t) + tx{t) (42) 



we need to choose 

m = ^ - M(i) + m (43) 

To test this scheme we numerically solved Eqs. ^ and 
()10|) for the same nonlinear system that was considered in 
Section fV Ai carrying out simulations with both a fixed 
two-photon detuning and a time-dependent two-photon 
detuning given by Eq. (|43p. The behavior of the atom 
laser linewidth over time is shown in Figure O It is clear 
that using a time-dependent detuning prevents the dele- 
terious chirp of the line center, and recovers the under- 
lying narrow spectrum of the laser output. 
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FIG. 9: Comparison of linewidth over time for a condensate 
with significant nonlinearities that becomes depleted. Solid 
line shows standard outcoupling; dashed line shows the effect 
of our chirp correction scheme. Parameters: A'o ~ 10^, i^t = 
150rads-\ a = 4 x 10"" m, fco = 3.2 x 10*^ m~\ 



C. Pumping 

A pumped atom laser operating at steady state would 
experience no chirp, as the energy of the lasing mode 
would be stable by definition. Lasers operating well 
over threshold also experience mode-selection effects that 
help provide a stable mode, although previous work sug- 

Pts that it may be nontrivial to operate in this regime 
■ IT^. Pumping can also induce gain- narrowing to com- 
bat the linewidth-broadening effects of quantum noise. 
A continuous pumping scheme for BECs has not been 
demonstrated, however, and consideration of quantum 
noise contributions to the linewidth of an atom laser re- 
quires a model that goes beyond the semi-classical ap- 
proximation. Examination of these quantum effects in 
zero-dimensional models has occurred in a variety of con- 
texts [lO, [U, [g^l , but no model has examined the com- 
petition between the multimode effects and the quantum 
noise in these devices. A multimode quantum model 
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to investigate the quantum noise contribution to the 
linewidth of an unpumped atom laser will be the sub- 
ject of a forthcoming paper. 

VI. CONCLUSIONS 

We have examined the linewidth of experimentally re- 
alistic non-pumped atom lasers in a variety of regimes. 
In strong outcoupling regimes, output spectra show a va- 
riety of undesirable features. If the condensate does not 
remain single mode due to non-Markovian effects caus- 
ing back action of the beam on the condensate, then the 
spectrum of the output beam will also be multimode, 
with peaks corresponding to the energy spacing of the 
excited modes in the condensate. We have shown that is 
possible to use Raman outcoupling and momentum con- 
servation rules to selectively filter out some of the excited 
modes of the condensate, resulting in a cleaner beam. 

In the limit as the coupling strength becomes weaker, 
we show that the linewidth of non-pumped atom lasers 
has two main limits. Weakly outcoupled, single-mode 
atom lasers with negligible nonlinearities will approach 



the Fourier limit, where the linewidth is given by the in- 
verse of the out-coupling rate. This outcoupling rate is 
a function of the atomic properties as well as the Rabi 
frequency of the change of state of the atoms. When the 
condensate has significant nonlinearities, as is the case 
is most current experimental atom lasers, depletion of 
the condensate causes the chemical potential to decrease 
over time, resulting in a downwards "chirp" in the cen- 
tral frequency of the output beam during the outcoupling 
process. This effect causes significant broadening of the 
linewidth of the beam. If high flux is required, and there- 
fore a significant fraction of the condensate must be out- 
coupled, we demonstrate that it is feasible to use a chirp 
compensation scheme such as sweeping the detuning of 
the state-changing outcoupling process. 
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